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Abstract 

The uniformly valid approximation to solutions of the radial Schrodinger equation with 
power-law potentials are obtained by means of the explicit summation of the leading con- 
stituent WKB series. 
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Similar questions arise when we use the perturbation theory or the variational method in 
order to solve the radial Schrodinger equation 

^ = ™. W , Q W ^ W + ^-*). (D 

How to find the unperturbed Hamiltonian or how to find the trial function for an arbitrary 
given potential? The universal answers are absent. In this sense both mentioned methods are 
not complete. At the same time the WKB approximation is directly determined by a given 
potential. The WKB approach deals with a logarithmic derivative 



Y(r) = dln ^ T \ = exp ^ Y{r')dr' 



which satisfies the nonlinear Riccati equation 

^ + <y H >-f>. (2) 

In this approach two independent solutions y ± (r) of the Riccati equation are represented by 
their asymptotic expansions 

F±(r) = h- 1 (±Q'/ 2 + f] h n Y±(r)j (3) 

in powers of Plank's constant fi. The usual WKB approximation contains a finite number of 
leading terms Y^[r) from the complete expansion Y^(r). The conventional WKB approximation 
is not valid at turning points. 

As it is well known, the WKB series is divergent. The direct summation of a divergent series 
does not exist. By summing one means finding a function to which this series is the asymptotic 
expansion There are several investigations on properties of the WKB terms 0, 9. In recent 
years many studies have been devoted to extracting some useful information about the exact 
eigenfunctions from the divergent WKB series (see for instance [Q] and references therein). In 
our previous work || we reconstructed the WKB series as a sum of new constituent (partial) 
asymptotic series. The explicit summation of the leading constituent series was performed. As 
a result the new quasi-uniform approximation to logarithmic derivatives was obtained. This 
approximation was derived formerly [|| by other means. Our approximation reproduces the 
known satisfactory approximation [|J near turning points. The proposed method was verified in 



the case of the one-dimensional Schrodinger equation without singularities Q. In the present 
work we apply our approach to the radial Schrodinger equation with the centrifugal term which 
is singular at the origin (r = 0). 

The second-order quasi-uniform approximation has the following form 

Y ap (r) = Y(r;t) = &i(r)j/i(a;t) + h(r)y 2 (a;t) (4) 

where 

1 Q(r) 1 Q> Q" 

^ = W*WrW^ hl{r) = W\oW Ur)= w- (5) 

Here we introduce notations for logarithmic derivatives 

yi(a;t) = In (Ai(a) + tBi(a)) (6) 
da 

of linear combinations of well studied Airy functions Ai(a) and Bi(a) |||. The value of an 
arbitrary parameter t determines a selected particular solution. We also use notation 

y 2 (a;t) = ^[-8a 2 ( yi (a;t)) 2 - 4a yi {a;t) + 8a 3 - 3] . (7) 

We consider power-law potentials 

V(r) = ar k , a > 0, k > 0. (8) 

When I = we deal with the usual one-dimensional problem for odd states which was solved 
formerly. When I 7^ we have a new problem. In this case there are two turning points r_ and 
r + (r + > r_) which satisfy equation Q(r) = 0. Turning points separate three regions. 

In the first region where r > r + we must describe only the decreasing solution of the 
Schrodinger equation. Therefore in this case we select the unique particular solution Y(r;0) 
of the Riccati equation. 

In the second region where r_ < r < r + we must describe oscillating solutions of the 
Schrodinger equation. Therefore in this case we select two particular solutions Y(r; ±i) of the 
Riccati equation. 

In the third region where r < r_ we must reproduce the known behaviour of the exact 
solution ^ ex (r) of the Schrodinger equation when r — > 0. In this case we have following simple 
formulas 

Q(r)-> * 2 * ( * a +1) , *es(r) - r l+ \ Y ex (r)^ l -±±. (9) 
At the same time we can derive relations 

a(r)^ ai = ^— -L) , h(r)^ l -, b 2 r) -» — . (10) 

V 4 / r r 

Thus, the approximate logarithmic derivative tends as follow 

Y{r . t) ^ c -M. (ii) 

r 

where 

c(Z, t) = - (2a iyi ( ai ;t) + 3y 2 (a /; t)) . (12) 

Naturally we must demand 

c{l,t) = l + l. (13) 
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It is the quadratic algebraic equation for determining the value of t. One of two solutions leads 
to an unphysical singularity. The second solution (nearest to zero ) ti gives dependence of t on 
I. We present two numerical examples t\ = —0.00393839, t 2 = —0.000523169. Hence, in third 
region we select the unique particular solution Y(r;fy) of the Riccati equation. 

Matching particular solutions at turning points we obtain the continuous approximate wave- 
function fy ap (r) which is represented by following formulas 

ty 1 (r) = Ncosfaexp^- j Y(r';ti)dr'^ (14) 



when r < r_ 




tf 2 (r) 

xcosl/ v ' ' ' v ' ' dr' -4>i\ (15) 



when r_ < r < r + and 



*s(r) = £(-!)» exp ( f r+ Y (r';+i)+Y(r>; i) \ ^ . , >, {/ , : ())f// , , (1(J) 
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when r > r + . Here 



ir+ 



7T 

<fo = - - arctan (17) 



and N is a normalization constant. 

We have the new quantization condition 

;-|i ^+')-^-) d ,. =r( „ + i )+ ,„ n=0jl , 2 ... (18) 

which determines the spectral value E sp (n, I) of energy implicitly. We denote the approximate 
wavefunctions with E = E sp (n, I) as ty ap {r\ n, I). 

So the approximate eigenfunctions are defined completely. But the question arises about the 
optimal approximate eigenvalues because the value E sp (n, I) is not unique choice. 

Since explicit expressions for wavefunctions are already obtained we are able to calculate 
expectation values 

E(n, I) =< tf ap (n, l)\H r \V ap (n, I) > (19) 

of the radial Hamiltonian 

where \^ ap (n, I) > is the vector in the Hilbert space which corresponds to the function ^ ap {r; n, I). 

In accordance with the eigenvalue problem H\ty > — E\^ >=0we construct the discrepancy 
vector 

\D(e;n,l) >= H r \^ ap {n,l) > -e\V ap (n,l) > (21) 

with an arbitrary parameter e while H r and \^ a p(n,l) > are given. It is natural to require 
that the discrepancy vector should not contain a component proportional to the approximate 
eigenvector. In other words we consider the orthogonality condition < ty ap (n,l)\D(e;n,l) >= 
as a criterion for selection of the optimal approximate eigenvalue. As a result we get just the 
expectation value E(n, I) while E sp (n, I) does not fulfil the above requirement. It should be also 
noted that the scalar product < D(e; n, l)\D{e; n, I) > is minimized at e = E(n,l). 
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Now we must verify our approximation numerically for the exactly solvable problems. We 
shall compare the normalized approximate wavefunctions \& ap (r;n, I) with the normalized exact 
wavefunctions *& ex (r;n, I). 

In order to estimate the closeness of two functions /i(r) and f%{r) we consider two corre- 
sponding vectors \f\ > and I/2 > in the Hilbert space. Then we construct a deviation vector 
I A/ >= |/i > — 1/2 > and a scalar product 

< A/|A/ >=< h\h > + < h\h >-< h\h >-< /2I/1 > • (22) 
Now we can define the relative deviation 

6f=1 - <h\h> + <h\h> (23) 

as the numerical estimation of the closeness of two functions. Note that Sf = when f\(r) = 
/ 2 (r). 

Thus, we get the following estimation 

5*(n,l) = 1- < * ex (n,l)\* ap (n,l) > (24) 

in the case of normalized real functions *$> ex (r;n,l) and *$> ap (r;n,l). The same numerical com- 
parison may be performed 

2<V ex (n,l)\% p (n,l)> 



6t>'{n,l) = 1 



(25) 



< K x (n,l)\% x (n,l) > + < % p (n,l)\% p (n,l) > 

for first derivatives ^'(r) = d^>{r)/dr. Naturally, we can define analogous estimations for higher 
derivatives. 

Now we compare two functions H r ^ ap (r;n,l) and E(n,l)^f ap (r;n,l). This comparison may 
be performed when we do not know exact solutions. As a result we get the relative discrepancy 



d(n, I) 



< V ap (n,l)\H?\y a Jn,l) > -E(n,l) 



< ^ ap (n,l)\H 2 \^ ap (n,l) > +E(nJ) 2 

which is directly connected with the initial Schrodinger equation under consideration. 
Finally, we characterize our approximation by the usual relative energy error 



(26) 



8E(n,l) 



E(n, I) 
E ex (n,l) 



1 



(27) 



where E ex (n, I) is the exact energy value. 

Firstly, we estimate our approximation in the case of the oscillator potential V(r) = a>2r 2 
for which the exact wavefunctions are well known ||. Table 1 shows that the quasi-uniform 
approximation gives fairly accurate results for all considered quantities in the oscillator case. 

Table 1. V(r) = a 2 r 2 . 



n, I 


5V(n,l) 


<5*'(n,0 


d(n, I) 


SE(n,l) 


0,0 
1,0 
2,0 


1.595 10" 5 
1.504 10~ 6 
4.271 10~ 7 


7.118 10" 5 
3.098 10~ 6 
6.391 10~ 7 


6.087 10~ 5 
7.312 10~ 7 
8.413 10~ 8 


5.635 10" 5 
1.884 10~ 6 
3.348 10~ 7 


0,1 

1,1 

2,1 


1.132 10~ 3 
1.337 10~ 3 
1.506 10- 3 


6.625 10~ 3 
2.898 10~ 3 
2.891 10~ 3 


9.475 10~ 3 
7.719 10~ 4 
5.863 10~ 4 


3.083 10~ 3 
1.011 10~ 3 
7.271 10" 4 


0,2 
1,2 
2,2 


5.052 10" 4 
4.173 10" 4 
4.635 10" 4 


4.877 10~ 3 
8.886 10" 4 
8.646 10" 4 


4.627 10~ 3 
1.153 10~ 4 
8.713 10~ 5 


1.492 10~ 3 
2.487 10" 4 
1.813 10" 4 
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The second example is the linear in r potential V(r) = a,\r . When I = the quasi-uniform 
approximation reproduces the exact result. When I ^ the exact wavefunctions are unknown 



but there are results of numerical solution for some energies 1C ] . We calculate values of d and 
5E for these states. Table 2 demonstrates validity of our approximation in the linear case. 



Table 2. V{r) 



ot\r 



n, I 


d(n, I) 


5E(n,l) 


0,1 


6.951 10" 


-3 


2.684 10" 3 


1,1 


5.357 10" 


-4 


8.230 10" 4 


0,2 


2.565 10" 


-3 


1.160 10~ 3 


1,2 


8.652 10~ 


-5 


2.238 10~ 4 



Hence, the performed reconstruction of the WKB series and subsequent explicit summation 
of the leading constituent (partial) series yield the satisfactory (qualitative and quantitative) 
description of wavefunctions in the case of the radial Schrodinger equation with power-law 
potentials. 
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